A Possibility of Search for New Physics at LHCb 
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It is interesting to search for new physics beyond the standard model at LHCb. We suggest that 
weak decays of doubly charmed baryon such as Scc(3520)''' , Sf,^ to charmless final states would be 
a possible signal for new physics. In this work, we consider two models, i.e. the unparticle and Z' 
as examples to study such possibilities. We also discuss the cases for Hj^j, which have not been 
observed yet, but one can expect to find them when LHCb begins running. Our numerical results 
show that these two models cannot result in sufficiently large decay widths, therefore if such modes 
are observed at LHCb, there must be a new physics other than the unparticle or Z' models. 

PACS numbers: 14.80.-j, 13.30.Eg 



I. INTRODUCTION 

LHC will begin its first run pretty soon, and besides 
searching for the long-expected Higgs boson, its main 
goal is to explore new physics beyond the SM. Many 
schemes have been proposed to reach the goal. Indeed, 
the LHCb detector, even though is not responsible for the 
Higgs hunting, will provide an ideal place to study heavy 
flavor physics and search for evidence of new physics. 
One can make careful measurements on rare decays of 
B-mesons, b-baryons, B-mixing and CP violation with a 
huge database available at LHCb, moreover, we are in- 
spired by the possibilities of discovering new physics. It 
would be beneficial to conjecture more possible processes 
which would signal existence of new physics. 

In 2002, the first event for doubly charmed baryon, 
S+ (3520), was observed by the SELEX Collaboration in 
the channel of S+ A+A'-7r+ [J, [3, [3]. S+ has the 
mass TO = 3519 ± 1 MeV and width F < 5 MeV. By 
studying an alternative channel of pD^K~ conducted 
later, the mass of the baryon-resonance was confirmed 
as TO = 3518 ± 3 MeV [1], which is consistent with that 
given in Ref. [l|. In the present theory, there definitely 
is no reason to exclude existence of which contains 
ecu valence quarks and as well {bhu) and SjJ'j^ {bhd), 
by the flavor-SU(3) symmetry. 

In this work, we propose that direct decays of Sec with 
charmless final states or with bottomless final states 
would be signals for new physics. By the quark-diagrams, 
one can easily notice that the main decay modes of Sec 
would be D+K{YP), A^K^, D+PK~ and KcK-it+. The 
later two modes are just the channels where the SELEX 
collaboration observed the baryon Sec. While the di- 
rect decays of Sec (or Sh;,) into charmless (bottomless) 
final states are suppressed in the standard model, so that 
would be sensitive to new physics beyond the SM. 

Since in Sec there are two identical charm quarks which 
can neither annihilate, nor exchange W-boson to convert 
into other quarks. In the SM, direct transition of See into 
charmless final states may realize via the double-penguin 
mechanism which is shown in Fig. 1 (a), the crossed 
box-diagram (Fig. 1 (b)) and a possible two-step process 



shown in Fig. 1 (c). The mechanism includes two pen- 
guin loops or a crossed box-diagram is very suppressed, 
so that cannot result in any observable effects and we can 
ignore them completely. If a non-zero rate is observed at 
LHCb, it should be a signal of new physics. Definitely 
the diagram of Fig. 1 (c) may cause a non-zero contri- 
bution and contaminate our situation for exploring new 
physics. If we consider the charmless decays of S++ or 
bottomless decays of S^j^, that diagram (Fig. 1 (c)) does 
not exist at all. Then, the first question is that can we 
distinguish such direct decays of Sec into charmless final 
states (or Sf,6 into bottomless final states) from the sec- 
ondary decays which result in charmless (or bottomless) 
products and are the regular modes in the framework of 
the SM. The answer is that the direct transitions are fa- 
vorably two-body decays, namely in the final states there 
are only two non-charmed hadrons by whose momenta 
one can re-construct the invariant mass spectra of See 
(or Sbb), whereas, in the regular modes with sequent de- 
cays, there are at least three hadrons in the final states. 

The second question is that is there any mechanism 
beyond the standard model available which can result 
in such direct decays? Below, we use two models to 
demonstrate how such direct decay modes are induced 
and estimate the widths accordingly. One of them is 
the unparticle scenario and another one is the SU (3) x 
SU{2)l X SU{2)u X U{1)b-l model where a new gauge 
boson Z' exists and mediates an interaction to turn the 
charm quark into a u-quark. Thus by exchange of an un- 
particle or Z' between the two charm quarks in 'Rtc^^'^ 



these direct 



(or between the two bottom quarks in ^^j, 
transitions occur. 

In this work, for simplicity, we only consider the inclu- 
sive decays of S++ (^^b) i'^to charmless (bottomless) final 
states. The advantage of only considering the inclusive 
processes is obvious that we do not need to worry about 
the hadronization of quarks into final states because 
such processes are fully governed by the non-perturbative 
QCD effects and brings up much uncertainty. 

Below, we will investigate the processes caused by ex- 
changing unparticles and Z' separately and then make a 
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FIG. 1: (a) The double-penguin diagram which can induce the 
decay of E.cd'S.bb) into non-charm (non-bottom) final states, 
(b) The crossed box diagram, (c) An emission where the effec- 
tive interaction would be non-local and for charmless decays 
of S'^^ , it does not exist. 



brief discussion on the possibility that new physics may 
result in observable phenomena at LHCb. 



II. THE INCLUSIVE DECAY OF DOUBLY 
CHARMED BARYON 

A. The unparticle scenario 

Before entering the concrete calculation, we briefly re- 
view the concerned knowledge on the unparticle physics 
0, which is needed in later derivation. The effective La- 
grangian describing the interaction of the unparticle with 
the SM quarks is 



A; , 



(1) 



where Ou and 0[J are the scalar and vector unparticle 
fields respectively, q and g' denote the SM quark fields. 
Generally, the dimensionless coefficients d^^^ is related to 
the flavor of the quark field. This interaction induces a 
FCNC and contributes to the processes of concern. 



For a scalar unparticle field, the propagator with mo- 
mentum "p and scale dimension du is 



with 



Ad.. ^ 



2sin(d;^7r) (— p2 - ie)2-<iM 

167r5/2 T{du + 1/2) 
(2^)2'i" T{du-l)n2duy 



(2) 



(3) 



with du the scale dimension. 

For the vector unparticle, the propagator reads 



d4xe'P-^(0|TO^O^(0)|0) 



(4) 



where the transverse condition d^O^ 



is required. 

In the unparticle physics, the inclusive decay of doubly 
charmed baryons into light quarks ccq — > uuq occurs at 
the tree level, and the transition is depicted in Fig. [51 
Here the exchanged agent between the two charm quarks 
can be either scalar or vector unparticle. 



FIG. 2: The inclusive transition of doubly charmed baryon in 
unparticle physics, where the double-dashed line denotes the 
scalar or vector unparticle in the unparticle model or Z' in 
the left-right model. 




Even though we only consider the inclusive processes 
where the quarks in the final states arc treated as on-shell 
free particles and the wavefunctions of the light hadrons 
in the final states are not needed, the binding effect of the 
initial baryons (Sj^"*" or 2^) which are composed of three 
valence quarks must be taken into account. Namely, 
when we calculate the hadronic matrix elements, we need 
to invoke concrete phcnomcnological models to carry out 
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the computations where the wave function of the initial 
baryon is needed. In this work, we adopt a simple non- 
relativistic model, i.e. the harmonic oscillator model 01 ■ 
This model has been widely and successfully employed in 
similar researches Ii,li,[ifl[Il[il[il[il. Thus one can 
trust that for heavy hadrons, such simple non-relativistic 
model can work well and the results are relatively reli- 
able, even though certain errors are not avoidable. Thus 
in this work the matrix elements of the effective opera- 
tors evaluated in terms of the harmonic oscillator wave 
function is believed to be a good approximation. Ac- 
cording to the references listed above, the errors in the 
estimate, especially as we only need the wavefunction of 
the the initial hadron, are expected to be less than 10%. 
By changing the input parameters and the model param- 
eters which are obtained by fitting other experiments, we 
scan the region of changes of the numerical results and 
find that the error range is indeed consistent with our 
expectation. 

In the harmonic oscillator model, the wave function of 
the doubly charmed baryon Sec (3520)+ is expressed as 

color. spill 

xV'H+(Pp,PA)6t(pi,s;)6t(p^,4)6t(P3,4)|0), 
which satisfies the normalization condition 

(*^+(P,s)|vI/^+(P',s')) 



(27r 



i3 ^ 



Ep 



■5\V -V')6{s- s'), 



where M is the normalization constant. Xsj^ ^^"^ 
note the spin-flavor and color parts of the wavefunction 
of doubly charmed baryon Scc(3520)+ respectively whose 
explicit expressions are 




2\c]c]di)^\c]cid])-\c[c]d]) 



In the harmonic oscillator model, the spatial wavefunc- 
tion 7/1^+ (Pp,Pa) reads as 



V'h+ (Pp,Pa) 



P^_^ 

2ap 2ax 



and the parameters ap and ax reflect the non- 
perturbative effects and will be given in later subsection. 

In the center of mass frame of Sec (3520)"'", the hadronic 
matrix elements Sfi is written as 

% = i2n)U\pi + P2 + P3 - M)T 

withr= (Ts + Ty). 

For exchanging scalar unparticle, Ts matrix element is 
written as 



Ts = Y. j d%d^Px{27r) 



3 ^P3 



spin 



x[m„(pi, si)7^(l - jsWip'i, s[) 

XU„(p2, S2)7i/(1 - l5)Uc{P2, s'2)] 



2sin(dw7i') (-p^ - ie^-du 
xA^^Z-^-f (Pp,Pa). (5) 



A, 



For the vector unparticle exchange, Ty is 



Tv ^Y. j d^Ppd'pxi^Ti) 

spin 



3 -^P3 



X K(pi, si)7p(l - 75)uc(pi, s'l) 

XU„(p2,S2)7i^(l - lb)Uc{p'2^ s'2)] 

f c^' ^2 z(~g^- + g^g7g2) 

A^"-i 2sm{dun) {^p^ - ie)^-du 

X-^V'b+ (Pp'Pa)- (6) 

Here Uq and Uq {q ~ c, u) denote the Dirac spinors 



' Eq + Vlq ( 1 



2mq V E^+m^ 



I En +m 



2mn 



X' 1, - 



(J ■ p 

Eq + niq 



and we can use the expression [15| 

jc^p 6m Ain\ sin di(TT\ 



|Cs"P 2m Am| sindj^Trj 



(7) 

(8) 

(9) 
(10) 



to simplify Ty and Tg. One needs to sum over all possible 
spin assignments for the Dirac spinors. 



with the definitions 



The Z' scenario 



/ / T->' ^ T->' „/ 

Pi P2 Pi + P2 - — P3 

Pp = Pa 



V2 

Pi+P2+P3. 



2 2mc-\-ma 
rrid 



The Left-Right models [l6| is also a natural extension 
of the electroweak model. It has been widely applied to 
the analysis on high energy processes. For example, re- 
cently He and Valencia [17[ employed this model with 
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certain modifications to explain the anomaly in ob- 
served at LEP [l^. Barger et al. studied Z' mediated 
flavor changing neutral currents in B- meson decays [l9l |. 
Bs — Bs mixing [23| and B Ktt puzzle [2l| . 

The gauge group of the model [l^ is SU{3) x SU{2) l x 
SU{2)ji X U{1)b-l where the four gauge couplings 
<73i 9L, QR and g correspond to the four sub-groups re- 
spectively. The vacuum expectation values of the three 
Higgs bosons break the symmetry. The symmetry break- 
ing patterns are depicted in literature. The introduction 
of a scalar field (j) causes Zq in the standard model to 
mix with a new gauge boson Zji, then Z, Z' arc the 
mass eigen-states. 

For the neutral sector the Larangian is 



Then we have the final expression as 



^ 2cos^6'w ^^^^^^ - 5A75)'7(cos^zZ^ - sin^zZ^) 
+^tan6'i^(iqL7^9L + i^umYum - '^dm^^dRi) 
(sin^zZ;, - cos^zZj, ) 

+ ^(tan0R + cot6iR)(sinezZ^ - cosCzZ;j 
{V^uVLjdml^dRj - VS,VR,^URa^uR,). (11) 



Here 9w is the electroweak mixing angle (tan^w = 
9r parameterizes the relative strength of the right- 
handed interaction(tan6'r{ = ^). is the Z — Z' mix- 
ing angle and V^j'^ are two unitary matrices that rotate 
the right-handed up-(down)-type quarks from the weak 
eigen-states to the mass eigen-states. Note that we use 
current notation for Pati-Salam model, and only third 
family couples to SU{2)r in this model. 

In the Z' model, inclusive decay of doubly charmed 
baryons into light quarks ccq — uuq occurs at tree level. 
The Feynman diagram (Fig. [5]) is the same as that for 
the unparticle scenario, but only the exchanged agent is 
replaced by Z' . 

In the center of mass frame of Sjl"^, we can obtain 



T = 



E 

spill 
X 



TTLp^ 



Uu{pi, si)7^(l - 75)iic(pi, s'l) 

XSu(P2,S2)7iv(l - j5)Uc{p'2, s'2)] 

ig^" r(7Ltan0w(tan0R -I- cot0R)cos^z 



(p2 + M'i) ^ 2 



(12) 



The authors of Ref. [l7|, 122, |23| suggested that cot6'R is 
large, so that tan^R can be ignored. They took approx- 



imations tan0wcot6'R^^ ~ 1 and cos^z ~ 1- Because 



E / dh,d^p,{2.f^ 

spin ^ 

x[u„(pi,si)7^(l -75)uc(Pi,s'i) 

XUuiP2, S2h„{l - 75)Wc(P2; 4)] 

x{-ignGFV2iVi^;,V^i,^f 
xAA^2+(Pp,Pa). (13) 



C. The expression of the decay width 

The inclusive decay rate would be obtained by inte- 
grating over the phase space which involves three free 
quarks and the procedure is standard [23], 

r(Scc(3520)+ -^uud) 

/•b2 /•2tt r-l IT'P 

(14) 

where ai, 02, &i and &2 are defined as respectively 



ai ~ 0, a2 
b 



2 2 M3+ 

1 = -^[cr(T + m+m_) - ^Jpf{T - m\){T -ml)], 



62 = ^[cr(T + m+m^) + \Jpf{T - m\){T - rn?.)]. 



Mz' is larger than 500 GeV [T] 



one has , , ^ 



1 



a = AI~+ — Pi, T = (T^ — Y (p5) , TO± = m-2 ± "m^. 

Here M-++, mi, m2 and m^ denote the masses of the 
doubly charmed baryon, up and down quarks respec- 
tively. In the following, for obtaining numerical results, 
we use the Monte Carlo method to carry out this inte- 
gral. For baryon S^, the expression is the same but only 
the mass of charm quark is replaced by that of bottom 
quark. 



III. NUMERICAL RESULTS 

Now we present our numerical results. 

Since only S++ has been measured, in the later cal- 
culation, we use its measured mass as input, and for 
^bb '^^ '^^^^ °^^y illustrate the dependence of its de- 
cay rate on the parameters. The input parameters in- 
clude: Gf = 1.166 x 10-5 GcV-2, ma = 1.60 GeV, 
mu = md = 0.3 GeV, m, = 0.45 GeV, m^ = 4.87 GeV. 

= 3.519 GeV. Up = 0.33 GeV2,aA = 0.25 GeV^ 
p] |25|. [26l . \2T\ . Here, the light quark mass refers to the 
constitute mass. 
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TABLE I: The decay widths of Hcc(3520)+ -» uud or S++ 
uuu and H^^ — > ddd{ssd) corresponding to du = 3/2 (in units 
of GeV). In the table, the second, third and fourth columns 
respectively correspond to the contributions from exchanging 
scalar unparticle, vector unparticle and both. 



FIG. 3: (a) and (b) respectively show the dependences of 
the decay widths of Hcc(3520)''' uud and HjJ'^ — > ddd re- 
spectively coming from the contributions of scalar unparticle, 
vector unparticle and both on du- 





scalar vector scalar -|-vector 


r[Hcc(3520)+ -> uud] 


4.57 X 10"^® 1.11 X 10"^^ 1.24 x 10"^^ 


r[E-, ^ ddd] 


5.85 X 10"^° 1.65 X 10"^^ 1.83 x 10"^^ 




5.21 X 10"^° 1.23 X 10"^^ 1.44 x 10"" 



A. The results in the Unparticle scenario 

For the unknown parameters An in the unparticle sce- 
nario, according to the general discussion, the energy 
scale may be at order of TeV, thus one can fix Au = 1 
TeV. We choose du = 3/2 in our calculation. 

In this work, we also calculate the inclusive decay 
width of doubly bottomed baryon S^^. The mass of ST^ 
is set as 10.09 GeV according to the estimate of Ref. [23], 
although there are no data available yet. 

The numerical results are provided in Table ID Fig. 
[3] illustrates the dependence of the decay widths of 
Scc(3520)"'" uud and ddd on du, the three lines 

(solid, dashed and dotted) correspond to the contribu- 
tions of scalar unparticle, vector unparticle and both on 
du- It is noted here "both" means that at present we 
cannot determine whether the unparticle is a scalar or 
vector and it is also possible that both scalar and vec- 
tor exist simultaneously. Thus wc assume both of scalar 
and vector contribute and they interfere constructively. 
Definitely, it is worth of further investigation. 



■D 

"i 

(0 

u 

V 

Q 



All (10""GeV) 

Scalar (10""GeV) 

Vector (10 "GeV) 



1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 



(a) 



All (10""GeV) 

■ - - Scalar (10""GeV) 

■ ■ ■ Vector (lO'^'GeV) 




The results for Z' exchange 



(b) 



The earlier studies indicate that the mass of Mz> 
should be larger than 500 GeVjll and V^^^V}^;^ is 
bound no more than 2.0 x 10~^ |22l |. in our calculation, 
we take their extreme values as Mz' = 500 GeV and 
^RtcVRtu = 2.0 x lO"'^, thus we would obtain the upper 
limit of the decay width. It is estimated with all the 
input parameters as 



^[S,,(3520)^ 



uud] = 7.66 X 10"^' GeV. 



(15) 



This is a too small numerical value compared with the 
width of S+ , therefore, it is hopeless to observe a non- 
zero branching ratio of S+ into charmless final states if 
only Z' is applied. 



Estimate the contribution from Standard Model 



As indicated above, in the framework of the SM, S+ 
can decay into two-body final states via Fig. [T] (c). It 
would be interesting to compare the SM contribution 
with that from the two models. Thus, we would roughly 
estimate the ratio of the contribution of Fig. [1] (c) to 
that of Fig. [5] (b) for the Z' model. It is easier to com- 
pare them because the structures of two diagrams and 
the relevant effective vertices are similar. 

By the order of magnitude estimation and with the 
SU(3) symmetry, the ratio of the amplitude of Fig. [1] (c) 
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{Tsm) to Fig. m (b) (T„„ or Tz') is 



SM 



SG f 



,2 V ^[^Bii^Rij) 



(16) 



where g is the momentum of unparticle or Z' (in the 
case of the SM, <^ ^^w) ^^^"^ "-^"^ be neglected in the 
propagator. Because the contribution from smaller 
is dominant, in the estimation, we set = 0.5GeV^. 
Since the whole case under consideration, may fall in 
the non-perturbative QCD region, as a rough estimate, 
we take a, = 1, V^^.V^^^ = 2 x lO"" and Gp = 
1.166 X 10-5GeV"^ we can get ^ « 66 (i.e.^ « 

^ z' '- z' 

4300 ). Then we can obtain the ratio of decay widths 
^V"^-^-^""- ~ 40. This ratio indicates that for S+ the 
contribution of the SM is smaller than that of the unpar- 
ticle scenario, but larger than that from the Z' model. 

However, for 2^^^, Fig[TJc) does not contribute at all, 
so that the decay of into charmless final states (or 
S^f, into bottomless final states) is more appropriate for 
exploring new physics than S+ . 

It is worth noticing that the estimate of the contribu- 
tion of the SM to the decay rate is very rough, thus what 
we can assure to ourselves is its order of magnitude. In- 
deed the magnitude contributed by the SM is very small 
and cannot produce sizable observational effects at all, 
even though it has a comparable order with that from 
the two sample models, the unparticle and Z' . In the 
future, if such mode were observed at LHCb, we can def- 
initely conclude that it is not caused by the SM, but new 
physics. 



IV. DISCUSSION AND CONCLUSION 

In this work, we propose to explore for new physics 
beyond the SM at LHCb by measuring direct decays of 
S+ (5^j, S|J{,) into charmless (bottomless) final 

states. Such decays can occur via the diagrams shown 
in Fig. 1 in the framework of the SM, but is much sup- 
pressed to be experimentally observed, therefore if a siz- 
able rate is measured, it would be a clear signal for new 
physics beyond the SM. We use two models as exam- 
ples, namely the unparticle and Z' models to calculate 
the decay rates, because both of them allow a transition 
of cc (bb) — > qq where q may be light quarks to occur 
at tree level. Thus one expects that these new models 
might result in non-zero observation. 

Indeed, our work is motivated by three factors, first 
the great machine LHC will run next year and a remark- 
able amount of data will be available, then secondly, the 
double-charmed baryon Sec which was observed by the 
SELEX collaboration provides us a possibility to probe 
new physics, and the last reason is that some models have 
been proposed and they may induce a flavor-changing 
neutral current, concretely the unparticle and Z' mod- 
els are employed in this work. Definitely none of the 
two models are confirmed by either theory or experiment 



yet, and they still need further theoretical investigations, 
but their framework is clear, so that we may use them 
as examples to demonstrate how new physics may cause 
such decay modes and indicate that a sizable observa- 
tional rate is a clear signature for new physics beyond 
the SM. Moreover, the double charmed baryon has only 
been observed by the SELEX collaboration, but not at 
B-factories. It seems peculiar at first glimpse, but care- 
ful studies indicate that it is quite reasonable due to the 
fragmentation process of heavy quarks. The authors of 
Ref.fl^ indicate that the meson Be cannot be seen at 
any e+e~ colliders because its production rate at such 
machines is too small, but by contraries, its production 
rate is greatly enhanced at hadron colliders. It was first 
observed at TEVATRON and its production rate at LHC 
would be much larger by several orders [2^. In analog, 
one can expect that such double-charmed baryons Sec or 
double-bottomed Sf,b can only be produced at LHC, but 
not at B-factories. 

The inclusive decays of doubly charmed baryons 
Scc(3520)+, S++ and S°j,, Sj^^ are explored in unpar- 
ticle and Z' scenarios. Our result indicates that the up- 
per limit of the inclusive decay width of S+T*" uuu 
is about 10~^^ GeV with du = 3/2. For inclusive decay 
^bb ^ ddd{ssd), the upper limit is at order of 10~^^ GeV. 
It is learnt that in the unparticle scenario, the contribu- 
tion from exchanging a vector unparticle is much larger 
than that from exchanging a scalar unparticle, as shown 
in Table U 

The parameters which we employ in the numerical 
computations are obtained by fitting other experimen- 
tal measurements, for example if the recently observed 
£)0 _ jjo j^g interpreted by the imparticle model, an 
upper boimd on the parameters in the model would be 
constrained. Indeed, all the present experimental data 
can only provide upper bounds on the model parameters 
no matter what new physics model under consideration 
is. 

So far it is hard to make an accurate estimate on the 
production rates of the heavy baryons which contain two 
heavy quarks at LHC yet, but one has reason to believe 
that the production rate would be roughly of the same 
order of the production rate of Be which was evaluated by 
some authors [2^ , or even smaller by a factor of less than 
10. The production rates indeed will be theoretically 
evaluated before or even after LHC begins running. 

In Ref. [2^, the authors estimate the number of Sec 
produced at LHCb as about 10^. Since the available en- 
ergy is much higher than the masses of Sec and Sf,;,, one 
has strong reason to believe that their production rates 
are comparable. Unfortunately our numerical results in- 
dicate that the unparticle and Z' scenarios cannot result 
in sizable rates for S++ uuu — > two hadrons and 
'^bb ^ ddd{ssd) —^ two hadrons which can be measured 
at LHCb and neither the SM. Even though the two sam- 
ple models and SM cannot cause sufficiently large rates, 
the channels still may stand for a possible place to search 
for new physics. If a sizable rate is observed at LHCb, it 
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would be a signal of new physics and the new physics is 
also not the unparticle and/or Z', but something else. 
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